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Abstract. In view of the role of reaction equations in physical problems, the 
authors derive the explicit solution of a fractional reaction equation of gen- 
eral character, that unifies and extends earlier results. Further, an alternative 
shorter method based on a result developed by the authors is given to derive 
the solution of a fractional diffusion equation. Fox functions and Mittag-Leffler 
functions are used for closed-form representations of the solutions of the respec- 
tive differential equations. 

1 Introduction 

Fractional reaction and diffusion equations involve fractional derivatives with re- 
spect to time and space and are studied to describe anomalous reaction and dif- 
fusion of dynamical systems with chaotic motion. Fractional reaction equation 
for Hamiltonian chaos is discussed by Zaslavsky (1994). Solutions and applica- 
tions of reaction equations are studied by Saichev and Zaslavsky (1997). Solu- 
tions of a fractional reaction equation is investigated by Haubold and Mathai 
(2000) for a simple production-destruction mechanism. This equation was gen- 
eralized by Saxena, Mathai, and Haubold (2002). In recent articles, Saxena, 
Mathai, and Haubold (2002, 2004a, 2004b) discussed the solution of a number 
of generalized fractional reaction equations. In the present article, we investigate 
the solution of a unified fractional reaction equation, which provides unification 
and extension of results on fractional reaction equations given earlier by Haubold 
and Mathai (2000) and Saxena, Mathai, and Haubold (2002, 2004a). We also 
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present the solution of a fractional integral equations discussed by Miller and 
Ross (1993). Further, an alternative proof of the solution of a fractional diffu- 
sion equation given earlier by Kochubei (1990) is investigated, which is based 
upon a result given by Saxena, Mathai, and Haubold (2006). Most of the re- 
sults are obtained in terms of generalized Mittag-Leffler functions in elegant and 
compact forms, which are also suitable for numerical computation. 

The paper is organized as follows. Section 2 provides the solution of a 
unified fractional reaction equation while Section 3 considers special cases of 
the equation. A shorter alternative method for the solution of a fractional 
diffusion equation discussed earlier by Kochubei (1990) is presented in Section 
4. A series representation and asymptotic expansion of the solution are given 
in Section 5. An H- function representation of a one-sided Levy stable density 
is also obtained. 



2 Fractional reaction equation 

In this Section, we present a method based on Laplace transform for deriving 
the solution of the unified fractional reaction equations. 

Theorem 1. If Re{vj) > 0,aj > 0,j £ N, and f(t) be a given function, defined 
on SR_|_, then the equation 

n 

N{t) - No fit) = -J2^j oOr^Nit), (1) 
is solvable and its particular solution is given by 

m = E ^ (..),,"'L,)i {n'°'-'"} 

1=0 ri + ...+r„_i=( Kft=l ) 

where the summation in (2) is taken over all nonnegative integers ri, . . . ,r„ 
such that ri -|- . . . -|- r„_i = I, and provided that the series and integral in (2) 
arc convergent. Here qD^ ,j € N are Riemann-Liouville fractional integrals, 
defined by 

oA""^/(i) = -pfr /*(* - ur-'f{u)du, Re{v) > 0, (3) 
i \^) Jo 

with oD°f{t) = f{t) (Oldham and Spanier, 1974; Miller and Ross, 1993), 
E^^^{z) is the generalized Mittag-Lefiler function, defined by Prabhakar (1971) 
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in terms of series representation as 

Proof. By the application of the convolution theorem of the Laplace transform 
(Erdelyi et al., 1953, p. 259) to (3), we find that 

L{oI?r/W;s} = L|^|L(/(i)), 

= s-\r{-s), (5) 

where /""(s) = e^''*' f{t)dt,s £ C,Re{s) > 0. Applying Laplace transform 
to (1) and using (5), it gives 

N-is) = - (6) 

If we employ the identity (Abramowitz and Stegun, 1968, p. 823) 
ix,+... + Xmy= E (rJ^^\r )i n^^'^' 

ri+...+r„=; ^=1 

where the summation is taken over all nonnegative integers, ri,...,r„, such 
that ri + . . . + r„ = Z, then for \ais~'^^ \ < 1, (7) transforms into the form 



=0 n+...+r„_i=, 



{n::;i(«M+i)'^''} 



ri>...r^_l>0 



(i + ais--o'+' ■ 

Taking the inverse Laplace transform of (8) by making use of the formula (Kil- 
bas, Saigo, and Saxena, 2004, eq. (12)) 

i-i {s-7(i _ as-^)-^; t} = t-^-^El^^iat^), (9) 

where Re{s) > \a\^^^'^'-'^\ Re{j) > 0, i?e(s) > 0, and applying the convolution 
theorem of the Laplace transform, the result (2) is established. 

RemEirk 1. The generalized Mittag-Leffler function defined by (4) is studied by 
Prabhakar (1971) and Kilbas, Saigo, and Saxena (2004). Recently this function 
is used in the theory of finite-size scaling of systems with strong anisotropy and 
long-range interaction by Chamati and Tonshev (2006). 
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3 Special cases 

Some special cases of Theorem 1 are of interest to be highlighted. If we set 
Vj = jv,aj = i^)c'''{j S N), we obtain 

Theorem 2. If Re{v) > 0, c > and f{x) G K+, then the equation 

N{t) - Nofit) = - f^(;?)c-Z?,— iV(i), (10) 

r=l 

is solvable and its solution has the form 

N{t) = Noj^ f{u)E:^,[-c''{t - ur]du, (11) 

where E" i (x) is the generalized Mittag-Leffler function defined by (4) and pro- 
vided that the integral (11) is convergent. 

When n = 1, we obtain the following result given by Hille and Tamarkin 
(1930). 

Corollary 2.1. Let Re{v) > 0, c > and let f{x) £ 3?+, then for the solution 
of the integral equation 

N{t) - Nofit) = -c^ oD^^Nit), (12) 

holds the following formula 

N{t) = Noj^ f{u)E,[-c-{t - ur]du, (13) 

where E^{z) is an entire function of order p = ^ and type a = 1, defined by 

oo ^ 

^-('^) = E Yi^^TTTy ^ ^' ^^^""^ > ^^^^ 

Note 1. The above result has also been given by the authors in a different form 
(Saxena, Mathai, and Haubold, 2004a, 2004b). 

If we set f{t) = t'^-^Ef,^^[-{cty], Theorem 2 yields 

Corollary 2.2. Let Re{y) > 0,i?e(7) > 0,c > 0, then for the solution of the 
reaction equation 

n 

N{t) - Not-^-'Et^^i-ictr] = - J2(rV oD^^^n G N (15) 

r=l 

holds the relation 

N{t) = 7Voi^-'<+"[-(ci)1,n e N. (16) 
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For f{t) = tP ^, Theorem 2 yields the fohowing result 

Corollary 2.3. If Re{p) > 0, Re{u) > 0, c > 0, then for the solution of the 
equation 

n 

N{t) - Not"-' = - J2irV oDr''N{t),r e N, (17) 
holds the relation 

N{t) = N^t>'-^V{p)Eiy{ctrir ^ N. (18) 

For n = 1, eq. (18) reduces to a result given by Saxena, Mathai, and Haubold 
(2002, p. 283, eq. (15)). When = a^s"^, for j = 1, . . . , u, we obtain 

Theorem 3. Let Re{i^) > 0,a > 0,t > 0,n > 1, \a"+'^s-'-'^+^>\ < 1, and f{x) 
be a given function defined on JR+, then the equation 

n 

Nit) - Nofit) = -J2a^ oD^'^^Nit), (19) 

r=l 

is solvable and its solution is given by 

mt) = 7Vo|^^V(«)i^(n+i).,.K+'(i-«)^"+'^'']c!« 

-aj\t - «)'^-ii?(„+i).,,K+i(t - u)^"+'>]du^ , (20) 
where E(^n+i)u,i'i^) is the generalized Mittag-Leffler function Ea,piz) defined as 

oo ^ 

Eoc,p iz) = E p(^^ ^ ^) , (a, ^ G C, Reia) > 0, ReiP) > 0) (21) 

and provided that the integral in (20) is convergent. 

If we take = ju, for j = 1, . . . n, then it is interesting to note that Theorem 
1 yields the following result given by (Miller and Ross, 1993) in a different form: 

Theorem 4. Let Re{v) > 0,aj > 0, and fix) be a given function defined on 
3ff+, |ais~''| < 1, then the fractional reaction equation 

n 

Nit) - Nofit) = - E «i oD;^-Nit), (22) 
is solvable and has the solution given by 

Nit) = No±{-iy (.,)■ Z ,), (n(«M+ir4 (23) 

X rf(u)it-u)'^:~>^''^'^'^^-'E^'+^„_ [-a,it - u)'^^]du, 

provided that the series and integral in (23) are convergent. 
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4 Fractional diffusion equation 

In this Section we present an alternative shorter method for deriving the solution 
of a diffusion equation discussed earlier by Kochubei (1990). 

Theorem 5. Consider the Cauchy problem 



with 



o£'"iV(x, t) = -c^ANix, t), (0 < a < 1; a; € < t < T), (24) 



N{x, t = 0) = S{x), X e lim 7V(x, t) = (25) 

a:— f ±C30 



of" is the regularized Caputo (1969) partial fractional derivative with respect 
to t, defined by 



QD^N{x,t) 



1 



r(i - a) 



d N{x,s)ds N{x,0) 



dt 



(t - xY 



and A is the Laplacian. The fundamental solution of the above Cauchy problem 
is given by 



N{x,t) = \x\-''tt-^H{ 



2,0 



^^u l(n/2,l),(l,l) 



(26) 



where h1'2{-) is the H-function (Mathai and Saxcna, 1978). 



Proof. Applying the Laplace transform with respect to t, using the result 
(Caputo, 1969) 

m— 1 

L{QD^N{x,t)) = s"N'-{x,s) - s°'-'-'^N^''\x,0), 

r=0 

m — I < a < m,, m S N, 
and Fourier transform with respect to x, gives 

s"iV'^*(A;,s) -s"-^ = -c''\k\'^N'-'{k,s), 

where the symbol " ~ " indicates the Laplace transform with respect to the time 
variable t and the symbol " * " the Fourier transform with respect to the space 
variable x. 

Solving for N'" , we have 



N~ {k,s) 



(27) 



By virtue of the following Fourier transform formula (Samko, Kilbas, and Marichev, 
1990, p. 538, eq. (27.1)) 
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(F,[\x\^'--y'K^r.-2)/2{a\x\)]){r) 



2tt 



n/2 



-^,(t G SR";n e N,a > 0), 
(28) 

where the multidimensional Fourier transform with respect to a; e SR" is defined 

by 

{F^N){T,t)= [ N{x,t)e'='''dx (t GW;t>0) (29) 
and ifi/(.) is the modified Bessel function of the second kind, yields 



iV(a;,s) = c-'^s"-'(27r)-2 



\X\C2 

s"2 



-2/2 



S 2 X 



(30) 



In order to invert the Laplace transform, we employ the following result 
given by the authors (Saxena, Mathai, and Haubold, 2006) 





(p,2<t) 
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(31) 



where K^{x) is the modified Bessel function of the second kind, Re{z'^) > 
0, Re{s) > 0. Thus we obtain the solution in the form 



N{x,t) 



H 



2,0 
1,2 



) 2 C 2 


"i" x 2^ 5 




(l-f-^,a) 


Ac"" 





By virtue of a result in Mathai and Saxena (1978), 



|(op,ap) 



,{ap+<7Aj,,Ap) 
\{bg+aBg,Bg) 



(32) 



(33) 



the power of the expression [{t "jxp} /Ac''] can be absorbed inside the H- 
function and consequently we obtain 



N{x,t) = \7rix\-" Hll^ 



t-"\x, 
Ac" 



(f,l),(l,l) 



Remark 1. If we employ the identity (Mathai and Saxena, 1978) 



Tjm.n 



A|(ap,Ap) 



Tjm.n 



|(ap,Ap/A) 



,A> 



(34) 



(35) 



the solution given by (32) can be expressed in the form 



1 



N{x,t) = ^k*a;r" hI'° 



(4c-) J 



(1,1) 

(f.^).(i.^). 



(36) 
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where a > 0. We also note that the above form of the solution is due to 
Schneider and Wyss (1989). There is one importance of our result (32) that it 
includes the Levy stable density in terms of the H- function as shown in (34). 
Similarly, using the identity (33) we arrive at 







2ct 


(t.^).(i,^). 



(37) 



where n is not an even integer. This form of the H- function is useful in deter- 
mining its expansion in powers of x. Due to importance of the solution, we also 
discuss its series representation and behavior. 



5 Series representation of the solution 

Using the sciries expansion for the H-function given in Mathai and Saxena (1978), 
it follows that 



tt2,0 
-"1,2 



X 



(1,1) 



J_ f r(t-f)r(i-^) 

2m J, r(l-s) "^'^^ ^^^> 



a 



^ r(i-f -0(-i)'x"(t+') ^ r(f - i-0(-i)'x"(i+') 

^ r(l - 3^-al){l)\ r(l-a-aO(/!) 



where n is not an even integer. 
Thus for n = 1, we find that 

1 °° A~ 
^(-'*) = ^E(-l)' r(l-a(/ + l)/2)(Z!) ' (^9) 

where A = ^ and the duplication formula for the gamma function is used. 

For n = 2, the H-function of (37) is singular and in this case, the result is 
explicitly given by Barkai (2001) in the form 

Nix,t)r^ ln[^-l]. (40) 

ttI (1 — ajt" x 

For n = 3, the series expansion is given by 

^(^' *) = 47ri3«/MV2 E r[l-a(l + Z/2)] ' ^^^^ 
From above, it readily follows that for n = 3 and a ^ 1 

N(x,t) ~ -, as X ^ oo. (42) 

X 
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It will not be out of place to mention that the one sided Levy stable density 
can be obtained from Laplace inversion formula (31) by virtue of the identity 



(43) 



and can be conveniently expressed in terms of the Laplace transform 



Jo 



/ e-"*$p(t)dt = e-"^ Re{u) > 0, Re{p) > 



(44) 



The result is 



^p{t) = -Hi:', 



i (1,1) 



{P > 0). 



(45) 



Note 2. This result is obtained earlier by Schneider and Wyss (1989) by fol- 
lowing a different procedure. Asymptotic behavior of $a(i) is also given by 
Schneider (1986). 

In conclusion, we mention that some of the results derived in this article 
may find some applications in problems associated with models of long-memory 
processes driven by Levy noise and other related problems, see the article by 
Anh, Heyde, and Leonenko (2002). 
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